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Abstract. With the recent accelerating rate of discoveries in the field of nearby young stars, the 
ability to identify new nearby young stars is as important as ever, and membership identification 
codes will continue to perform a vital role in scientific research. In the process of creating a 
new moving group membership identification code - LocAting Constituent mEmbers In Nearby 
Groups (LACEwINC) - we have compiled a few pointers relevant to astronomers trying to use 
codes like LACEwING to locate young stars. 
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1. Introduction 

One of the most important methods used in characterizing young stars, and the only 
method that can establish membership in specific nearby young moving groups, is their 
kinematics. Kinematics exploits the property that young stars (at ages less than a billion 
years) are still tracing the space motion of the gas from which they formed. Several 
membership identification codes are available publically: BANYAN (Malo et al. 2013), 
BANYAN II (Gagne et al. 2014); a convergence code (Rodriguez et al. 2013); and (soon), 
LACEwING (Riedel et al. in prepjj- 

These codes use the six basic kinematic elements (RA, DEC, tt, ^ira, f^DEC, and RV) 
to predict probabilities of membership in nearby young moving groups, and can handle 
incomplete data. It is, however, important to know how to use and interpret the results 
of LACEwING (and codes like it) to obtain the most accurate results. 


2. Overview of LACEwING 

The LACEwING code computes membership probabilities in 10 moving groups (e 
Chameleon, TW Hydra, (3 Pic, Octans, Tucana-Horologium, Columba, Argus, AB Do- 
radus, Hercules-Lyra, and Ursa Major) and 4 open clusters [rj Chameleon, The Pleiades, 
Coma Berenices, and The Hyades). The groups are represented as freely oriented triaxial 
ellipsoids, which are computed from new lists of bona-fide (high confidence) members 
taken from existing lists (e.g. Malo et al. 2013), and filtered using epicyclic tracebacks 
to limit membership to only the stars that were plausibly within the group’s boundaries 
at the time of formation. 

The LACEwING code matches stars to groups by computing up to four metrics of 
membership, matching the proper motion, distance, radial velocity, and spatial location 
of the star to those estimated for a group member at the same RA and DEC. Those 

t See also https://github.com/ariedel/lacewmg 
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metrics are combined into a single numerical goodness-of-fit value. The goodness-of- 
fit values are transformed into membership probabilities by creating a large simulation 
(using 32 million stars; Riedel et al. submitted) of the Solar Neighborhood, which is then 
binned by goodness-of-fit value relative to a group. This is shown for /3 Pic and Coma 
Ber in Fig.[l] where all seven different possible combinations of data (RA and DEC are 
assumed to be known) are given as different colors. Curves are fit to those distributions 
to provide membership probabilities. 



Goodness of Fit 



Goodness of Fit 


Figure 1. Histograms of matches to the simulated P Pictoris Moving Group (left) and Coma 
Berenices open cluster (right), showing the percentage of stars in each bin that are genuine 
members of the group/cluster. Curves fit to these histograms (dotted lines) are used to translate 
the LACEwING code’s goodness-of-fit parameters to actual membership probabilities. 


3. Lessons from LACEwING 

3.1. Groups Are Not Equally Easy to Identify 

Some groups (like the Coma Berenices open cluster) have very distinctive UVW space 
velocities and space positions (Fig.[2|), which translate into uniquely identifying proper 
motion and radial velocity estimates for any given star. 

As an example, the large simulation of stars in Fig.[T] shows us that for the case in 
which all kinematic elements are available (RA, DEC, plus proper motion, distance, and 
rv), a goodness-of-fit value of 0.75 relative to (3 Pic yields a roughly 50% probability that 
the object is a member of /3 Pic. The stars matched to Coma Ber show that, for the 
equivalent case using all available data, a goodness-of-fit of 0.75 yields a nearly 100% 
probability that the object is a member of Coma Ber. 

3.2. Additional Kinematic Data Vastly Improves Membership Probabilities 

Fig. [T] shows the seven different data availability scenarios LACEwING expects: RA and 
DEC plus all combinations of proper motion, radial velocity, and distance. For (3 Pic, 
the proper motion method (red dashed line) demonstrates that even in the best possible 
case of a goodness-of-fit = 0, the probability of membership is only 35%, i.e. 65% of 
those perfect matches to /3 Pic are actually (simulated) members of other moving groups 
and the generic field population. Adding a radial velocity (pink triangles) increases the 
best-case probability to 80%/20% contamination. 

Another good reason to obtain more kinematic information is to reduce the false pos¬ 
itive contamination (Fig. [3]) of a survey. With only proper motion and a fixed cutoff of 
20% probability, over 80% of the objects in a survey will be contaminants (to be weeded 
out by some other technique, such as spectroscopy). Having complete data for the star 
of interest reduces the contamination rate to just above 40%. 
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Figure 2. 2D projections of the freely oriented triaxial ellipsoids that constitute LACEwING’s 
kinematic models (in UVW motion and XYZ spatial position) of the 14 nearby young mov¬ 
ing groups and open clusters. As shown in these figures, some groups have overlapping UVW 
space velocities or XYZ space positions. This can make it difficult to determine to which group 
individual stars belong. 



Figure 3. In this figure we see the cost of lowering the membership probability threshold: the 
false positive contamination of the sample increases. The best way to reduce the false positive 
contamination is to obtain more of the 6 kinematic elements for the sample of stars, but even in 
this case, additional methods (spectroscopic confirmation, for example) will be needed to weed 
out the sample. 

3.3. Not All Young Stars Belong to a Known Moving Group 

Virtually all lithium is primordial and easily fused by objects over roughly 60 Jupiter 
masses. Therefore, a low mass object with detectable lithium must either be a brown 
dwarf, or a young star. Riedel et al. (submitted) assembled a sample of roughly 500 
lithium-detected nearby stars and used those stars as a test sample for LACEwING, 
reasoning that these should be the best nearby young moving group members. However, 
nearly half the sample did not match any of the known nearby young moving groups 
LACEwING tests for. The same lithium-detected sample was run through other moving 
group codes - BANYAN I and II, and the Convergence code - with similar results. Some 
lithium-detected stars were not sorted into moving groups (Fig.|4]). 

There are three potential reasons for this: 1.) There are additional nearby young stellar 
groups not currently known or accounted for; 2.) there is a young held population not 
part of any previously identihed groups; 3.) our current understanding of the kinematics 
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Figure 4. 2MASS 0447-7412 has a significant lithium equivalent width, 50±20 mA(Murphy & 
Lawson 2015), a proper motion, a radial velocity, and no distance. If we generate a range of 
reasonable distances for the star, it could spatially be within the Octans (cyan) moving group, 
but it cannot possibly share the space velocity of Octans, and cannot be a member (as previously 
noted by Murphy) 

of moving groups is incomplete. The fact that multiple moving group codes (each with a 
different algorithm and kinematic model of the nearby young moving groups) could not 
identify large numbers of stars as members of nearby young moving groups suggests that 
the explanation probably involves some combination of the first two points. Ultimately, 
the message is clear: Young stars do not need to be members of a known nearby young 
moving group. “None of the above” can be a correct answer. 


4. Conclusions 

As we enter the third decade of nearby young star research, there is much more to be 
done from the standpoint of membership statistics. Moving group identification codes 
such as LACEwING will continue to provide essential tools in studies of stellar and 
planetary formation and evolution, but care must still be taken in interpreting their 
results. 
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Discussion 

Shkolnik: What about chemical abundance tagging? 

Riedel: LACEwING does only one thing: it uses kinematic and spatial data to evaluate 
potential memberships in moving groups. It should be used alongside spectroscopic meth¬ 
ods to confirm that the star(s) in question is/are actually young, and have appropriate 
ages for the group LACEwING suggests is the most probable host. 

Zuckerman: Why isn’t Carina-Near in the code? 

Riedel: There wasn’t enough in the literature about it. The same goes for Cas-Tau, 32 
Ori, and Alessi 13. But they can be added fairly easily - new ellipsoids, new simulation. 












